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Abstract 

We demonstrate that a transition from decelerated to accelerated cosmic expansion arises as a 
pure interaction phenomenon if pressureless dark matter is coupled to holographic dark energy 
whose infrared cutoff scale is set by the Hubble length. In a spatially flat universe the ratio of the 
energy densities of both components remains constant through this transition, while it is subject 
to slow variations for non-zero spatial curvature. The coincidence problem is dynamized and 
reformulated in terms of the interaction rate. An early matter era is recovered since for negligible 
interaction at high redshifts the dark energy itself behaves as matter. A simple model for this 
dynamics is shown to fit the SN la data. The constant background energy density ratio simplifies 
the perturbation analysis which is characterized by non-adiabatic features. 
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I. INTRODUCTION 



Nowadays an overwhelming, direct and indirect, observational evidence supports the idea 
that the Universe is currently undergoing a phase of accelerated expansion. The most recent 
and precise confirmation of this rather unexpected feature (see, however [l, 2]) was provided 
by the 3rd year data of the WMAP mission 31 . Likewise, data from high-redshift supernovae 

n n n 

type la the cosmic microwave background radiation [5j, the large scale structure j6|, the 



integrated Sachs- Wolfe effect p], and weak lensing |8|], endorse it to the point that the 
discussion has now shifted to when the acceleration began and which might be the agent 
behind it. For recent reviews see j^. 

According to our understanding on the basis of Einstein's gravity, more than 70% of the 
cosmic substratum, dubbed dark energy, must be endowed with a high negative pressure. 
Less than 30% is found to be pressureless matter. Most of the latter is in the form of cold 
dark matter; only about 5% corresponds to "normal" baryonic matter. Currently, we neither 
know what the dark matter is made of nor do we understand the nature of dark energy. 

Lacking a fundamental theory, most investigations in the field are phenomenological and 
rely on the assumption that these two unknown substances evolve independently, i.e., their 
energies are assumed to obey separate conservation laws. In particular, this implies that 
the dark matter energy density varies as a~^, where a is the scale factor of the Robertson- 
Walker metric. The behavior of the density of the dark energy is then entirely governed by 
its equation of state, the determination of which is a major subject of current observational 
cosmology. One should realize, however, that any coupling in the dark sector will change 
this situation. A coupling will modify the evolution history of the Universe. As one of 
the consequences, the energy density of the (interacting) dark matter will no longer evolve 
as a^^. Ignoring a potentially existing interaction between dark matter and dark energy 
from the start, may result in a misled interpretation of the data concerning the dark energy 
equation of state. Das et al. [l^ and Amendola et al. [llf have shown that a measured 
phantom equation of state may be mimicked by an interaction, while the bare equation of 
state may well be of a non-phantom type. Further, models showing interaction fare well when 
contrasted with data from the cosmic microwave background 12| and matter distribution 
at large scales 1^. It can therefore be argued that the possibility of dark energy to be in 
interaction with dark matter must be taken seriously. 
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On the other hand, there exist hmits for the strength of this interaction for various 
configurations [l^. It is characteristic for these approaches that they admit a non- interacting 
hmit. If the interaction is switched off, they continue to represent models of a mixture of dark 
matter and dark energy. The best known examples are models with a decaying cosmological 
"constant": if the decay is switched off, they reduce to the ACDM model. An interaction 
here leads to (possibly important) corrections of the non-interacting configuration. E.g., for 
a given (not necessarily constant) negative equation of state parameter of the dark energy, 
the interaction manifests itself in third order in the redshift in the luminosity distance of 

n 

supernovae type la [15]. It does not, however, (directly) influence the leading orders. 

The qualitatively new aspect of the present approach is the circumstance that an interac- 
tion is crucial already in leading order. In this context the accelerated expansion itself is a 
consequence of the interaction. The coupling does not just lead to higher order corrections 
of a non-interacting reference model which by itself provides an accelerated expansion of 
the Universe. The non-interacting limit of the present approach reproduces an Einstein-de 
Sitter universe, i.e., there is no accelerated expansion if the interaction is switched off. This 
non-interacting limit is supposed to characterize our Universe at high redshifts. 

The data also suggest that the present Universe is nearly spatially fiat Many 
studies neglect the spatial curvature term and focus solely on the spatially fiat case, thereby 
taking for granted that the spatial curvature term necessarily leads to trifling corrections. 
Several of the not so many works that retain that term, for the sake of generality, even 
seem to justify that it is of minor importance. On the other hand, thanks to the increasing 
observational precision also higher order corrections may become within reach in the not too 
distant future. The authors of p/7i] have demonstrated that the spatial curvature enters the 
luminosity distance of SNIa supernovae in third order in redshift (see also jisl). It is also 
known that there are degeneracies between the curvature and the dark energy equation of 
state in the corresponding parameter space [191]. The possible relevance of spatial curvature 
'or the lowest multipoles in the cosmic microwave background radiation was discussed in 
20| . These examples indicate that apart from general theoretical grounds, the observational 
situation will require the inclusion of spatial curvature, even if its contribution is small, at 
some level of precision. Further, as demonstrated by Ichikawa et al. |2l|, depending on the 
parametrization of the equation of state parameter the present value of the spatial curvature 
can be as large as 0.2. Therefore, the additional dynamics provided by the curvature, should 
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certainly not be dismissed too quickly since it contains information which is needed to further 
restrict the cosmological parameter space. In this connection, Clarkson et al. [22] have shown 
that by excluding the spatial curvature when interpreting the empirical data one can incur in 
gross mistakes when reconstructing the equation of state of dark energy. We mention that an 
effective spatial curvature term also app ears as the result of an averaging procedure within 



the "macrosopic gravity" approach 
on cosmological observations see, e.g 
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for a further discussion of the effects of averaging 
24b. 

In this paper we consider pressureless dark matter in interaction with an unknown com- 
ponent which is supposed to describe dark energy. We neither specify the dark energy 
equation of state nor the interaction rate from the beginning. These two (generally time 
dependent) parameters will influence the ratio of the energy densities of dark matter and 
dark energy. The behavior of this ratio is crucial for the "conventional" form of the "co- 
incidence problem", namely: "why are the matter and dark energy densities of precisely 
the same order today" . In principle, matter and dark energy redshift at different rates. We 
show, that there exists a preferred class of dark energy models for which the dynamics of the 
energy density ratio is entirely determined by the spatial curvature. For vanishing curvature 
the energy density ratio remains constant. These models are singled out by a dependence 
px oc H^, where px is the dark energy density and H = a/a is the Hubble parameter. 
Exactly this dependence is characteristic for a certain type of dark energy models, inspired 
by the holographic principle 25|, 1261] . 

Holographic dark energy models must specify an infrared cutoff length scale [s^]. The 
choice of this scale is presently a matter of debate. The most obvious choice, the Hubble 



length, seemed to be incompatible with an accelerated expansion of the Universe 



see, 



however. 
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30|). This is why, starting with the work of Li 31i, many researchers have 



adopted the future event horizon as the cutoff scale as this choice allows for a sufficiently 
negative equation of state parameter and hence an accelerated expansion -see, e.g. [32|]. 

In a previous paper [ssl, we showed that a cutoff set by the Hubble length may well be 
compatible with an accelerated expansion provided that the dark energy and dark matter do 
not evolve separately but interact, also non-gravitationally, with each other. In this setting, 
a negative equation of state parameter that gives rise to accelerated expansion arises as 
a direct consequence of the interaction. Here we put this feature in a broader context 
and demonstrate that a constant or slowly varying (as the consequence of a non-vanishing 
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spatial curvature) energy density ratio is compatible with a transition from decelerated to 
accelerated expansion under the condition of a growing interaction parameter. We show that 
in holographic dark energy models with a Hubble length cutoff a transition from decelerated 
to accelerated expansion is realized as a pure interaction effect. At high redshifts and for 
negligible interaction the dark energy equation of state approaches the equation of state 
for matter, such that a preceding matter era is naturally recovered. In this context the 
coincidence problem can be rephrased as follows: "why is the interaction rate between dark 
matter and dark energy of the order of the Hubble rate precisely at the present epoch?" 
This reformulation allows us to address the coincidence problem as part of the interaction 
dynamics in the dark sector. 

Different interaction rates will imply a different perturbation dynamics. We shall demon- 
strate that non-adiabatic features appear as a characteristic signature for a large class of 
interacting models. These effects may be used to discriminate between different models of 
the cosmic medium which share the same background dynamics. 

The outline of this paper is as follows. Section [III provides the general formalism for an 
interacting two-component fluid, where one of the fluids is pressureless. Then it focuses 
on the case for which the ratio of the energy densities of both fluids is constant or slowly 
varying. This will single out models for which the energy density of the second component 
is proportional to the square of the Hubble parameter. A realization of this dependence is 
provided by certain holographic dark energy models. The basic properties of these models are 
recalled in Section Unl Section ITVl discusses in detail matter in interaction with holographic 
dark energy. It provides the conditions under which an interaction driven transition from 
decelerated to accelerated expansion is feasible. It also comments on the possibility of 
a weak time dependence of the saturation parameter of the holographic bound, which is 
usually assumed constant. Section |V] introduces a simple interacting model and compares 
it with the ACDM model. The perturbation dynamics for a fluctuating interaction rate is 
discussed in|VT]with special emphasis on non-adiabatic aspects on large perturbation scales. 
Finally, section IVIII summarizes our findings. 
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II. INTERACTING COSMOLOGICAL FLUIDS 



The field equations for a spatially homogeneous and isotropic universe are the Friedmann 
equation 

3if2 = 87rGp-34, (A; = +1,0,-1), (1) 

and 

H =-A'kG {p + p) + ^ . (2) 
A combination of both equations yields 

a relation that will be useful later on. The deceleration parameter, q = — a/(aif^), can be 
written as 

'40 + ^)0+3^) ■ w 



The total energy density p is supposed to split into p = pm + Pxi where pm is the energy 
density of pressureless dark matter. Under this assumption the total pressure equals the 
dark energy pressure, p = px- We further assume that both components do not conserve 
separately but interact with each other in such a manner that the balance equations take 
the form 

Pm + ^HpM = Q and px + 3iJ(l + w)px = -Q , (5) 

where w = px/px is the equation of state parameter of the dark energy, and the function 
Q > measures the strength of the interaction. Models featuring an interaction matter- 



dark energy were introduced by Wetterich 34 1 (see also (35l]) and first used alongside the 
holographic dark energy by Horvat 29|] . Nowadays there is a growing body of literature on 
the subject -see, e.g. [36| and references therein. Although the assumption of a coupling be- 
tween both components implies the introduction of an additional phenomenological function 
Q, a description that admits interactions is certainly more general than otherwise. Further, 
there is no known symmetry that would suppress such interaction and arguments in favor 
of interacting models have been put forward recently |37|. 

The quantity of interest for analyzing the coincidence problem is the ratio r = Pm/px, 
which, upon using Eqs. ([5]), can be written as 
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r = (1 + r) 



3Hw 



1 + r 



+ r 



(6) 



Here we have introduced the quantity T = Q/px which characterizes the rate by which px 
changes as a result of the interaction. On the other hand, combining the balance equation 
for Px in ([5]) with Eq. ([3]), we obtain 



Px 



H 



3Hw 



1+r 



+ r 



+ 



1 + 3- 



w 



1 + r 



(7) 



Comparing now Eqs. ([7j) and it follows that the dynamics of the aforesaid ratio is 
governed by 



:i+r) 



k 



px 



1 + 3. 



w 



1 + r 



(8) 



It is this formula which deserves attention with regard to the coincidence problem. Clearly, 
the lower |r|, the less acute this problem seems to be. Inspection of Eq. ([8]) shows that the 
case Px oc is singled out for it leads to 



k I w 

1+r -TTT 1 + 3 



1 + r 



(9) 



At this point, it is expedient to introduce the dimensionless quantities 



n 



8tt G pm 



M 



Stt G Px 



3H^ ' ""^ 3H^ ' 
Thus, Friedmann's equation ([1]) can be cast as 

nM + nx + ^k = i- 

By using it together with Eq. (jl]), the expression iQ takes the form 

f = -2H-^q . 



(10) 



(11) 



(12) 



Via the curvature the evolution of r is directly linked to the deceleration parameter. For 
k = +1 and g > (decelerated expansion) r augments, and for g < (accelerated expansion) 



r diminishes. For k = —1 the behavior is just the opposite. In any case, only small variations 
of r are possible nowadays if H^) ^ 1 at present 0, [3]. While a slow time dependence 
of the energy density ratio is desirable from the point of view of the coincidence problem, 
it remains to be clarified whether dark energy models with px oc if ^ are able to account 
for a present phase of accelerated expansion as well as for a transition to the latter from an 
earlier matter dominated phase. 

A survey of the, by now, rather ample body of literature on dark energy candidates reveals 
that models of this type have indeed been introduced previously, namely, in the context of 
ideas which are rooted in the holographic principle. In the following section we briefly recall, 
how a dependence px oc emerges in these holographic dark energy models. 



III. HOLOGRAPHIC DARK ENERGY 



The basic ideas of the holographic principle were introduced by 't Hooft |25|] and Susskind 



261]. The development of interest for our purpose was put forward by Cohen et al. [2 
followed by Hsu [28] and Li 31|, who considered specific holographic dark energy models. 

The essential point in establishing the holographic idea is the way of counting the degrees 
of freedom of a physical system. Consider a three-dimensional lattice of spin-like degrees 
of freedom and assume that the distance between every two neighboring sites is some small 
length I which is of the order of the Planck length, ^pl. Each spin can be in one of two 
states. In a region of volume the number of quantum states will be = 2", with 
n = (L/i)^ the number of sites in the volume, whence the entropy, given by the logarithm of 
A^, will be oc (L/i)^ In 2. Identifying (in Planck units) with the ultraviolet cutoff A (the 
corresponding energy is L^A'^), the maximum entropy varies as S* ~ A^, i.e., proportional 
to the volume of the system. Based on considerations on black hole thermodynamics (bear 
in mind that the Bekenstein-Hawking entropy is Sbh = ^/ (4^p;), where A is the area of the 



black hole horizon), Bekenstein 38|] argued that the maximum entropy for a box of volume 
should be proportional to its surface rather than to its volume. In keeping with this, 't 
Hooft conjectured that all phenomena within a volume should be described by a set of 
degrees of freedom located at the surface which bounds this volume with approximately one 
binary degree of freedom per Planck's area. 



Inspired by these ideas, Cohen et al. 271] demonstrated that an effective field theory that 
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saturates the inequality (the Bekenstein bound) 

L^K^ <Sbh^L^MIi {Mil = {StiG)-^) , (13) 

necessarily includes states for which the Schwarzschild radius Rg is larger than the box size, 
i.e., Rg > L. Namely, for sufficiently high temperatures (in Planck units T ^ L^^ but 
T < A) the thermal energy of the system is ~ L^T^ while its entropy is S" ^ L^T^. When 
Eq. is saturated (by setting T = A in there) it follows that T ^ {MlJLf/^. Now, 

the Schwarzschild radius Rg is related to the energy by Rg ~ E/Mp^ (recall that an object 
with Schwarzschild radius Rg corresponds to a (Newtonian) mass Rg/{2G)). Consequently, 
Rg ^ L, i.e., the Schwarzschild radius is indeed larger than the system size. 

A stronger constraint which excludes states for which the Schwarzschild radius exceeds 
the size L is 

L^A^<Mj>iL. (14) 

The expression on the right-hand side of Eq. (|T4l) corresponds to the energy of a black hole 
of size L. So, this constraint ensures that the energy L'^A^ in a box of the size L does not 



3l|. 



exceed the energy of a black hole of the same size 

While the Bekenstein bound (fT3l) implies a scaling L oc A~'^, the inequality f|T^ cor- 
responds to a behavior L oc A~^. Furthermore, since saturation of f|T^ means A'^ ~ 
{MlJL^fl\ one finds that S = = {Ml^ L^f'^ = S\ 



■i/A 
BH- 



By saturating 



px it follows [271, 



;he inequality (031) and identifying A'^ with the holographic energy density 

3 



3c^ 

^^^^ = ' ^^^^ 

where the factor 3 was introduced for convenience and is a dimensionless quantity which 

is usually assumed constant. It is obvious, that for a choice L = H~^, i.e., for a cutoff set 

by the Hubble radius, the energy density (ITSl) is characterized by exactly that dependence 

oc which was singled out by Eqs. ([8]) and (Q. In the following we shall identify f|T5|) with 

the energy density px appearing in Eq. ([5]). 

IV. HOLOGRAPHIC DARK ENERGY IN INTERACTION WITH MATTER 

For holographic dark energy models with a cutoff scale H^^ the time dependence of 
the energy density ratio is determined by Eq. ([9]). For a spatially flat universe the ratio 



r remains constant. Possible changes of r due to a non- vanishing spatial curvature term 
are necessarily small. At first glance, models with a constant (or almost constant) energy 
density ratio seem unable to account for a transition from decelerated to accelerated 
expansion. And indeed, holographic dark energy models with the Hubble scale as cutoff 
length were considered to be ruled out since their equation of state parameter does not 
seem to allow negative values which are required for an accelerated expansion 28|, |31 |. 
However, as shown previously, the pre sence of an interaction between dark energy and dark 
matter may change this situation [33]. This becomes obvious if Eq. ([3]) with p = wpx is 
solved for the equation of state parameter w. The result is 

r ilk 



Assuming small, the interaction rate F essentially determines the equation of state 
parameter as soon as the dimensionless ratio T/H, which here and in the following is the 
relevant interaction parameter, becomes of order one. For negligible interaction, \T\/H ^ 1, 
one has \w\ <C 1. 

The main point of interest now is the behavior of the deceleration parameter, i.e., whether 
this kind of models admits a transition from decelerated to accelerated expansion. The 
condition for accelerated expansion (g < 0) 

w + —<--, (17) 
3H 3 ^ ^ 



readily follows from the above equations. Introduction of the expression (1161) for w into the 
latter leads to 

T ^ 1-nx _ ^ F^^ ^^g^ 

where we have used Friedmann's equation flTT]) . On the other hand, since holographic energy 
is not compatible with phantom fields [39;], one must impose w > — 1 to be consistent with 
the holographic idea. From Eq. (fT6|) we find that this condition amounts to 

^<3 (l-f]x)-^^fc. (19) 
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Therefore, the consistency condition simply reads, 

30x > 1 - fifc . (20) 
10 



The inequality (fT9i) says that the ratio on its left hand side is smaller than a quantity of 
order one. Clearly, this is compatible with an initial condition |r|/(3i/) ^ 1. The latter 
implies that w must be close to zero deep in the matter dominated era (cf. Eq. fll6p ). An 
initial condition |r| ^ i7 is equivalent to an initially negligible interaction. This means 
that at high redshifts (but well after matter-radiation equality), the dark energy component 
is almost non-interacting and its equation of state is close to that of pressureless matter. 
On the other hand, the condition f|T8|) for accelerated expansion entails that \T\/H must 
be larger than a lower bound of order one, as well. Both inequalities are consistent with 
condition ( l20l) . This supports our suggestion that, depending on the interaction rate, a 
transition from decelerated to accelerated expansion is indeed feasible. 

Starting from Eq. ([H]), the dynamics of the energy density ratio may alternatively be 
written as 

[L_ J . (21) 

For = we recover the stationary case r = constant, irrespective of the value of T/H. If 
additionally T/H = constant, it follows from Eq. ( |T6l) that w is also a constant. In such 



a case, there is no transition from decelerated to accelerated expansion [33|. However, if 
r/// is allowed to grow, then the transition may well occur, even though the ratio r remains 
constant or almost constant (see Eq. fl2^ below). 

The spatially curved cases generate an interesting additional dynamics. The transition 
condition T/H = r (cf. Eq. ( ITSi) ) coincides with the condition for r = for A; = ±1 in ( !2TI) . 
Since at early times one has T/ H <^ r, the ratio r in (12T]) increases for k = +1, while it 
decreases for k = —1. It reaches a maximum (minimum) for k = 1 [k = —1) at T/H = r 
and may go down (up) afterward. 

The obvious requirement for a sufficient growth of the interaction parameter T/H is that 
(T/Hr)' > in the phase of decelerated expansion or, by (|2T|) . that 



T/H 1-n 



X 



1 r' 



(22) 



But according to Eq. f|T9l) the ratio F/if is bounded from above. For a growth of T/H 
from a very small value to a value of order one, a transition from decelerated to accelerated 
expansion is feasible. We mention again, that any change of r arises solely as a curvature 
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effect. 

A different way to understand the dynamics of the ratio r, which also may be seen as a 
consistency check, is to start from Eq. (1111) and to reahze that Qx = constant. Then (recall 
that r = Qm/^x), 

'- — -IT,- (2^' 

The first term on the right hand side of fl23l) is constant. It is just the curvature term that 
makes r vary. Since in accelerating universes decreases r will approach a constant value 
in the long time limit. The ratio r grows if the curvature term enlarges the right hand side 
{k = 1 and g > 0, or A; = —1 and q < 0) and it goes down if the curvature term diminishes 
the right hand side {k = 1 and q < 0, or k = —1 and g > 0). This is essentially the behavior 
discussed beneath Eq. f|T2l) . Clearly, differentiation of Eq. fl23l) alongside the introduction 
of the deceleration parameter q consistently reproduces Eq. (fT2|l . 
Combining Eqs. (fT^ and fl2T]) we get 

77-'- = -2^'- (2^1 

This equation encodes the central message of our paper. It covers all the cases k = 0, ±1 
and explicitly shows that the cosmological dynamics crucially depends on the difference of 
two ratios: the ratio T/H and the ratio r of the energy densities. The difference of these 
two ratios is directly proportional to the (negative) deceleration parameter. Accelerated 
expansion requires T/H > r, decelerated expansion demands T / H < r. It is obvious, 
that even a constant or a slowly varying ratio r may be compatible with a sign change of q 
provided T/H evolves accordingly. Even with a constant r a present accelerated expansion is 
compatible with an earlier matter dominated period with decelerated expansion. The point 
is that for negligible interaction (at sufficiently high redshifts) the dark energy component 
behaves as non-relativistic matter. In this phase r is not really an important parameter 
since it describes the ratio of two components with the same equation of state. It was 
this property that apparently ruled out a (non-interacting) holographic dark energy model 



with an infrared cutoff set by the Hubble scale 28|, l3l|]- Here, this unwanted (in the non- 
interacting model) feature is advantageous since, thanks to it, a matter dominated phase 
during which structure formation can occur is naturally recovered. It is only because of the 
gradually increased interaction that the equations of state begin to differ from each other. 
Then r becomes important. At the first glance the circumstance that r is constant or only 
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slowly varying seems to imply, that the coincidence problem gets significantly alleviated. 
But in fact, it has been shifted from the problem to explain a constant (or slowly varying) 
ratio r towards the problem to explain an interaction rate which has to be of the order of the 
Hubble rate just at the present time. We argue that an alternative, dynamical formulation 
of the coincidence problem maybe useful, in particular, if a specific interacting model of 
the type presented here will lead to potentially observable differences from, say the ACDM 
model or generalized Chaplygin gas models (see below). 

One might argue against a constant -or nearly constant- ratio r by saying that, as is well 
known, at the time of primeval nucleosynthesis the dark energy should not contribute more 
than 5 per cent to the total energy density if the standard big bang scenario for the build 



up of light elements is to hold 



411 ]. However, this criticism does not apply to our case since. 



as said above, the equation of state of dark energy, w, was not close to —1 at early epochs 
but close to that of dust. 

By starting from the expression (fT5l) with L = H-^ for the dark energy we have assumed 
the parameter to be constant throughout our considerations. However, there does not 
seem to exist any compelling reason for the inequality ( |T^ to be saturated or for the degree 
of saturation to be a constant once and for all. In principle, a weak time dependence 
(0 < (c^) /c^ <ti H) of should be admitted here. Such an additional degree of freedom 
will modify the dynamics discussed so far. In particular, instead of described by Eq. fl2Tl) . 
the ratio r will evolve according to 



-Hr 



1-Qx 



1 r 

rH 



- 1 



+ 



(25) 



and the equation of state parameter changes from (fT6|) to 



w 



1-Qx 



3H 3 3H c2 



(26) 



As a consequence, in the conditions (fT7|) and (ITSl) and in all relations after them, the rate F 

(c2)' 

is to be replaced by F -|- . A varying will induce small variations in r, additionally to 
those which are due to the spatial curvature and can support the transition from decelerated 

n □ 

to accelerated expansion [33[ (see also I4(|). Thus, at this transition, we have 
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which reduces to the previous resuh r(g = 0) = for (? = constant. 

Before closing this section, it is noteworthy that a negative equation of state parameter 
w can be obtained from a negative curvature term {k = —1) alone (cf. Eq. f|T6|) ). i.e., even 

for r = 0, 



1 1 
w = — - 



3 a2i72 1 _ 



(28) 



However, the condition ( |T8l) for accelerated expansion, in this case l/{aH)'^ > 1 — Qx, can 
never be satisfied as it contradicts the constraint equation ffTTl) . Consequently, only models 
with a non-vanishing interaction between dark energy and dark matter can be potential 
candidates for a satisfactory cosmological dynamics based on the holographic idea with the 
infrared cutoff set by the Hubble function. 

V. A SIMPLE MODEL 

In our approach, the unknown nature of dark energy is mapped on a (so far) unspecified 
interaction in the dark sector quantified by the ratio T/H. In this section we build a simple 
model that explicitly exhibits the general features discussed in the preceding section. 

With the help of the equation of state parameter f|T6l) the balances expressions (|5]) may 
be written as 



Px 



1 

1 - - 



r Q 



X 



3 l-Qx 3H 1-n 



X 



P^ = -3H 

Pm 



1 - 



r VL 



X 



M 



(29) 



For k = 0, we have 1 — Qx = and both rates coincide. Moreover, in this case r stays 
constant. Integration of fl29l 2) under this condition from some initial period (subscript "z") 
onward, results in 



Pm 

PMi 







r daT 


1 exp 


r J 





(30) 
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The exponential describes the deviation from an Einstein - de Sitter universe. For the 
special case of a constant F/if, the expression (J30l) reduces to 



Pm 



( 



(31) 



pMi 



a 



Consequently, the Hubble parameter depends on the scale factor as H oc a r/(rH)] 
Hence, 



The exponent on the right hand side of the last expression is larger than unity for F/ (rH) > 1 
which reproduces the previous condition ( ITSl) for accelerated expansion. The limiting case 
r/{rH) = 3 corresponds to a constant energy density, i.e., to an exponential growth of the 
scale factor. 

As mentioned above, a constant ratio T/H can describe either a phase of decelerated 
expansion (V/{rH) < 1), or a phase of accelerated expansion {T/{rH) > 1), but never a 
transition between the two. To stage a transition from decelerated to accelerated expansion 
T/{rH) must increase which, for a constant ratio r, means that T/H must increase. Since 
there is no guidance of what a microscopic interaction model could be, we shall resort to a 
phenomenological approach, assuming a growth of the ratio T/H with a power of the scale 
factor a. Notice that a working microscopic interaction model in the present context, where 
it is exclusively the interaction that drives the accelerated expansion, would be equivalent to 
explaining the nature of dark energy. The fact that we obtain a transition from decelerated 
to accelerated expansion under the condition of a constant energy density ratio makes the 
coincidence problem appear in a different way. In our approach it remains to explain why 
the interaction rate is of the order of the Hubble rate just at the present time. We argue, 
that this reformulation of the coincidence problem may be advantageous since it offers a 
potential dynamical solution, albeit hypothetical at the present state of knowledge. 

In the following we assume a growth according to 



2 



a OC t3-r/('-ff) . 



(32) 





(33) 
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with constant, positive-definite parameters a and /3. The condition on /3 ensures, that 
the phantom divide cannot be crossed. Furthermore, we restrict our considerations to the 
past evolution of the universe. The ansatz fl33l) implies a continuous growth of the ratio 
T/{rH) with a maximum value at the present time. We do not speculate about a possible 
influence of the interaction on the future dynamics of the Universe. But it is obvious that 
if, in the future, the interaction becomes less effective, an evolution back to a fresh phase of 
decelerated expansion is feasible. 

Inserting the ansatz (133!) into Eq. (130!) we get 




The exponential factor diminishes the decreases of the energy density with cosmic expan- 
sion. Replacing the initial quantities by the corresponding present day values (subscript 
"0") provides us with 



PM = pMo \ —J exp 



.a \\ao 



(35) 



For k = and = constant, the energy density of the dark energy component shows 
exactly the same dependence on the scale factor. 



'aoY 
Px = Pxo [ —J exp 



3P_ 
a 



1 



(36) 



Thus, the Hubble rate is given by 



h = 



H 



a ) 



■i/2 



exp 



3^ 
2^ 



(37) 



Comparing Eq. (1331) with ( l24l) for the spatially fiat case, the deceleration parameter is 
found to be 



(38) 
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It is expedient to compare the dimensionless Hubble rate, h, given by Eq. fl37j) . with the 
corresponding quantity of the ACDM model, 



H 



ACDM 



PA 



pA + Pa/0 



_^ Pmo /^y 
Pa V a / 



1/2 



(39) 



Introducing the redshift parameter z hj ao/a = 1 + z, we have 



H 



ACDM 



PA 



Pa + Pmo 
Pmo 



2 Pa + Pmo 



l + ^{l + z)' 

PA 

z + 0{z^), 



1/2 



(40) 



where for reasons explained below we have retained only the linear term in z in the second 
line of last equation. From (1371) it follows that 




FIG. 1: Graphs of the dimensionless ratio h vs redshift. Except for the middle graph -which 
corresponds to the ACDM model-, from top to bottom the a values are: 2.0, 1.5, 1.2, and 1.0. 
Here, we have taken /? = 3/4. 
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H 



(1 + zf^^ exp 

-no 



3/3 



2a V V 1 + ^ 



= l + l{l-P)z + 0{z'), (41) 

for our interacting model. It is remarkable that to linear order in z, the expression ( HTi) does 
not depend on a. 

We may now require that our model (14T!) coincides with the ACDM model (140|) up to 
linear order in z. This allows us to fix the parameter (3 in terms of the observationally well 
established ratio of the parameters pmo and p\ of the ACDM model. The result is 

(3 = . (42) 

Pa + Pmo 

For the observed Pmq/pa ~ 1/3 we have P ~ 3/4. Under this condition, the present value 
of the deceleration parameter ( l38i) . (?o = | ~ f/^; coincides with the corresponding value 
Qacdmo = 1/2 — [3pa/2(pa + Pa/o)] ~ —0.6 of the ACDM model as well. However, the 
evolution of q towards this value depends on a which gives rise to differences in any order 
beyond the linear one. For the value of the transition redshift Zacc, that follows from q = 
in ( 138|) . we obtain Zacc = (3/9)"'^^" — 1. 

With the identification fH2]) and Pmo/pa ~ 1/3 we find Zacc ~ 1-2 for a = 1 and Zacc ~ 0.5 
for a = 2. The corresponding 2;-value for the ACDM model is approximately 0.8. 

Figure 1 contrasts the Hubble parameter prediction of our interaction model for different 
values of a with the ACDM model. Clearly, the free parameter a, obviously of the order of 
one, can be used to adjust h{z) according to the observational situation. 

Figure 2 presents the best fit of our model and the best fit of the ACDM model to the 
Gold SNIa data set of Riess et al. (fourth reference in ^), and the SNLS data set of Astier 
et al. (sixth reference in ^]). In plotting the graphs the distance modulus, p = 5 logdL + 25, 
was employed. In this expression = {1 + z) H~^{z') dz' is the luminosity distance in 
megaparsecs. As can be seen, both best fits largely overlap one another. In all, current SNIa 
data are unable to discriminate between the popular ACDM and our interaction model. 
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FIG. 2: Distance modulus vs rcdshift for the best fit model, 0,x = 0.75, Hq = 67.35 km/s/Mpc, 
a = 1.5 (solid line), and the ACDM model, 9,a = 0.71, and Hq = 67.35 km/s/Mpc (dashed line). 

VI. PERTURBATION DYNAMICS 



In a flat universe the dark energy equation of state parameter is given by 
w — — {1 + r) Since px — tt-, the total equation of state of the cosmic medium is 



P 
P 



3Hr 



(43) 



The corresponding adiabatic sound speed parameter becomes 



p 
P 



1 - 



1 iT/Hy\p 



3jH T/H / p ' 



7 = 1 + 



P 



(44) 



However, the interaction introduces non-adiabatic features into the perturbation dynamics. 
Deviations from the adiabatic behavior are suitably characterized by the quantity p — jp, 
which in our case takes the form 
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. p . 
p- -p 
p 



p 



V/H T/H p r 



(45) 



where a hat indicates perturbation of the corresponding quantity. The contribution on the 
right-hand side of (H5|) accounts both for the non-adiabaticity of the X component and for 
a contribution which results from the two component nature of the system. In terms of the 
components this expression is equivalent to 



P P \ Px J P^ Px 



D 



M 



(46) 



where 



D = 



P 



-3H 



P 



P + P P 

and the energy density perturbations for the components are 



Dm = -3H 



Pm 
Pm 



and Dx = —3H 



Px_ 
px 



Analogously, the pressure perturbations are described by 



P = 



P 



-3H- and P 



X 



-3H 



Px 



(47) 



(48) 



(49) 



P + P P Px 

The circumstance that r = in these class of models results in a particularly simple structure 
for Dx — Dm- Just from the definition of these quantities we find by direct calculation 



Dx-D 



M 



r 

7r 



(50) 



This difference is directly related to the fluctuation of the energy density ratio. The dynamic 
equation for Dx — Dm (see below) then determines the time dependence of f. 

A homogenous and isotropic background universe with scalar metric perturbations and 
vanishing anisotropic pressure can be characterized by the line element (longitudinal gauge, 



cf. 



[1 + 2^) dt^ + (1 - 2^) 5„^dx°da;^ 



(51) 



The perturbation dynamics is most conveniently described in terms of the gauge-invariant 



variable 
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P 



(52) 



which represents curvature perturbations on hypersurfaces of constant energy density. Cor- 
responding quantities for the components are 



Ca^-V'-^— (a = x,m). 

Pa 

On large perturbation scales the variable C, obeys the equation (cf. 



44 



45 



C 



-H\P-^D 
P 



The equation for Cx — Cm = \ {Dx — Dm) is 
(Cx-Ca/)' = 



Px 

pM Px 



p 



(53) 



(54) 



(55) 



The gauge invariant quantity Cn describes the perturbation of the interaction term. It is 
defined in analogy to the other perturbation quantities, 



n 
n 



n 



(56) 



Inserting in (|55l) the expression ( |50ll and 

i^/H) r 1 [T/H)- (p r \1 

V/H r(l + r) 37if T/H \p 1 + rJ ' ^ ^ 

which follows from the equation of state ( l43l) . we obtain 

f = =^ f = const (58) 

on large perturbation scales, i.e. scales for which spatial gradient terms may be neglected. 
This means, also at the perturbative level, f is not a dynamical degree of freedom on these 
scales. It is remarkable that this property holds for any interaction rate. While, of course, 
the trivial case f = is included here, it is expedient to notice that any non-vanishing f, 
although constant, gives rise to non-adiabatic effects. 

The simplest case to discuss this dynamics more specifically is F = const. Since H oc 
p^/^, this corresponds to an effective background equation of state p oc — p^^^ which is 



Px . Tp 
Px - —px = -WTT 
Px ^Hr 
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characteristic for a special generalized Chaplygin gas [47|. Its energy density is 

P = pi 



1 



(59) 



Again, for F = the Einstein de Sitter universe is reproduced. The general relations 
and fj45|) specify to 



P 



Ip 
2~p 



and 



. P . 
P- -P 
P 



1 



(0 = 3H) 



(60) 



(61) 



respectively. Only if spatial gradient terms in the perturbation dynamics are neglected, the 
first two terms in the bracket on the right-hand side of Eq. fl^Tl) cancel each other approx- 
imately. Therefore one rnay expect that non-adiabatic effects will be most important on 



small perturbation scales 48|, |49|. But even for - ^ 2^ a non- vanishing perturbation f gen 



erates a non-adiabatic dynamics of the system. The corresponding large scale perturbation 
dynamics can be solved analytically. The equation for ( in this case reduces to 



r f 

67 



(62) 



With 7 



1 — we may write 



where C = f • Since (cf. (EH])) 

H = H, 



C' = - 



r 



Tfl 



2 r a 3Hr — T 



(63) 



1 



r 



ai\3/2- 
a 



(64) 



equation (l63l) is readily solved. The result is 

. . rf 1 



3 r 3if,r — T 



3/2 



(65) 



The non-adiabatic contribution grows with a^^^. This will influence the time dependence 



of the gravitational potential and hence the integrated Sachs- Wolfe effect (cf 53|, |5J] for 



similar studies). Recall that the ACDM model is characterized by a constant value of ( 



On small scales we expect a modification of the adiabatic sound speed 



48 



49j. Notice that 



the quantity ? in fl60l) is negative. As a consequence, there occur small scale instabilities 
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within an adiabatic perturbation analysis of this generahzed Chaplygin gas. This apparently 
unrealistic feature has been used as an argument to discard this type of models altogether 



50l | . On the other hand, it has been suggested that non-adiabatic effects might modify this 



picture 



51 



52 



53| . That this is indeed the case, has been demonstrated recently by an 



explicit numerical analysis for the galaxy power spectrum |4^] . The present consideration in 
the context of interacting holographic models offer a systematic procedure to include such 
effects. We plan to perform a more quantitative non-adiabatic perturbation analysis in a 
subsequent work. 



VII. DISCUSSION 



We have established a holographic dark energy model in which a negative equation of state 
parameter as well as the transition from decelerated to accelerated expansion arise from pure 
interaction phenomena. A remarkable property of the model is the fact that, in a spatially 
flat universe, the ratio of the energy densities of dark matter and dark energy remains 
constant during the transition. The non-interacting limit of this model is an Einstein-de 
Sitter universe which is supposed to describe our cosmos at high redshifts. At flrst sight, 
the fact that r was never large may seem at variance with the conventional scenario of 
cosmic structure formation as one may think that at early times the amount of dark matter 
may have been insufficient to produce gravitational potential wells deep enough to lead 
to the condensation of the galaxies. However, this is not so; an earlier matter dominated 
phase is naturally recovered since for negligible interaction at high redshifts the equation of 
state of the dark energy is similar to that of non-relativist ic matter. In the context of our 
approach the coincidence problem is dynamized and takes the form: "why is the interaction 
rate that drives the accelerated expansion of the order of the Hubble rate precisely at the 
present epoch?" Based on the assumption that the relevant interaction parameter grows 
with a power law in the scale factor (Eq. fl33l) ). we have worked out a speciflc model for the 
transition from decelerated to accelerated expansion under the condition of a constant ratio 
of the energy densities of dark matter and dark energy. A preliminary analysis shows that 
this model flts the SNIa data not less well than the ACDM model. 

We have also shown that taking into account the spatial curvature term gives rise to an 
additional dynamics which implies a small (compared with the Hubble rate) change of the 
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energy density ratio. 

Furthermore, we discussed to what extent a slowly varying saturation parameter of the 
holographic bound may modify the cosmological dynamics. Besides, our approach can cope 



with a later transition to a new decelerated phase of expansion [55|] -something incompatible 
with holographic models whose infrared cutoff is set by the radius of the future event horizon. 

The interaction between dark matter and dark energy introduces non-adiabatic features 
into the perturbation theory. For the special case of a constant decay rate of the dark energy, 
in the background equivalent to a generalized Chaplygin gas, we find, Eq. (l65l) . that the 
large scale curvature perturbations on hypersurfaces of constant density (which are constant 
in the adiabatic case, in particular for the ACDM model) vary with the power 3/2 of the 
scale factor. 
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